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Abstract 

The connection between QCD, a nonlocal Nambu-Jona-Lasinio type model and the Landau 
gauge gluon propagator is explored. This two point function is parameterized by a functional 
form which is compatible with Dyson-Schwinger and lattice QCD results. Demanding the nonlocal 
model to reproduce the experimental values for the pion mass, the pion decay constant, T n ^^ and 
the light quark condensate we conclude that low energy physics does not distinguish between the 
so-called decoupling and scaling solutions of the Dyson-Schwinger equations. This result means 
that, provided that the model parameters are chosen appropriately, one is free to chose any of the 
above scenarios. Furthermore, the nonlocal Nambu-Jona-Lasinio quark model considered here is 
chiral invariant and satisfies the GMOR relation at the 1% level of precision. 
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The usual approach to Quantum Chromodynamics (QCD) defines the Green's function 
generation function a la Faddeev-Popov (FP) This approach allows the investigation of 
the perturbative solution of QCD. However, in what concerns nonperturbative phenomena, 
like confinement and chiral symmetry breaking, the FP trick is ill defined as it assumes a 
unique solution of the gauge fixing condition on each gauge orbit. Gribov demonstrated 
for the Coulomb and Landau gauges, that in each set of gauge related configurations there 
are more than one configuration, the Gribov copies, satisfying the gauge fixing condition. 
The study of Gribov was confirmed by many lattice simulations j^j, which also investigate 
how the copies modify the Green's functions. Furthermore, there is a no-go theorem 4] 
stating that the problem of the copies is common to all local gauge fixing conditions. 

Gribov tried to solve the problem of the copies by changing the functional integration 
space, replacing the hyperplane of transverse configurations, T = {A : d ■ A = 0}, by 
the first Gribov region Q: the subset of V with a positive defined Faddeev-Popov operator. 
This change of integration space modifies the tree level gluon propagator in the infrared 
region and, according to Gribov, at zero momentum the propagator vanishes. However, the 
problem of building a nonperturbative generating functional for QCD was not solved, as fl 
is not free of Gribov copies. 

In {5), |(J the restriction of the integration space to f2 was implemented via a modifica- 
tion of the QCD action. The Gribov-Zwanziger (GZ) action includes extra bosonic and 
fermionic fields, breaks softly BRST js, 7] and is perturbatively renormalizable. Moreover, 
the GZ action also predicts a vanishing zero momentum gluon propagator [§]. This result 



suggested a gluonic confinement criteria as discussed in [5 



a. 



Further extensions of the 



GZ action, including possible condensates, were investigated in [?| leading to the so-called 
refined Gribov-Zwanziger (RGZ) action. If the GZ action predicts a -D(O) = 0, according to 
the RGZ action D(0) ^ [10|. 

In the recent years there has been a considerable effort to compute the gluon and ghost 
propagators from first principles, namely using Dyson-Schwinger equations (DSE) and the 
lattice formulation of QCD. Considerable attention was given to the infrared propagators 
and to the value of D(0). 

DSE computations start from the Faddeev-Popov quantization procedure. The solutions 



follow into two categories: the scaling solution 
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12| which has D(0) = and is compatible 



with the Gribov-Zwanziger confinement criterium; the decoupling solution 



13Ml5| which 



predicts a -D(O) ^ 0, related to a dynamical generated gluon mass M(p 2 ). 



Wilson action 4D lattice QCD simulations for the 



g 



been performed for SU(2) and SU(3) gauge groups [16J422]. The gluon propagator has been 
computed for volumes as large as (27 fm) 4 for SU(2) 2l| (a ~ 0.21 fm), and (16 fm) 4 for 



uon and ghost propagators have 



SU(3) 22j (a ~ 0.17 fm). The raw data shows a decreasing D(0) with the lattice volume 
but no infrared suppression of Dip 2 ) has been observed. Furthermore, naive extrapolations 
of -D(O) to the infinite volume 20j do not provide a clear answer about the infinite volume 
value for D(0). 

In 2l|, inequalities relating the zero momentum gluon propagator with an average ab- 



solute value color-magnetization, which allow for a scaling analysis with the volume, were 
derived. SU(2) simulations give a finite and nonvanishing D(0) in the limit V — » +oo. 
A similar analysis for SU(3) 23j, although using smaller volumes, give -D(O) = in the same 
limit. The reasons for this difference are still not understood. 



In 



24j the authors considered ratios of lattice propagators, which, hopefully, suppress 



the finite volume effects. The analysis of the ratios point towards a vanishing -D(O) in the 
infinite volume limit. 

In conclusion, the present situation gives ambiguous results for the value of D(0). This 
question is of great theoretical relevance as it validates, or not, possible gluon confinement 
mechanisms. Any other way of estimating -D(O) independently is welcome. In this work we 
link -D(O) with low energy physics, i.e. the physics associated with pions. In order to connect 
the infrared gluon propagator with low energy phenomenology, an effective low energy chiral 
quark model of the Nambu-Jona-Lasinio type is built. 

In the QCD Lagrangian, the interaction between quarks and gluons reads 



(i) 



In the generating functional, expanding the exponential term containing C^ A to second 
order in g and integrating the gluon fields, the theory becomes an effective nonlocal fermionic 
theory. For one flavor QCD, after color and spin Fierz transformations, assuming that vector 
and axial vector currents play a secondary role in the dynamics, and taking into account 
only the color singlet bilocal currents J(x, y) = ip (x)^ l (y), Js(x, y) = ^(x) 75^(3/) (i stands 
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FIG. 1. Form factors as a function of p. The figure includes typical form factors used in previous 
studies - see, for example, 2g] and references there in. 



for color indices) one gets the effective action 



S[ip, ip]= d xd y\ ip(y) 8(y - x) (ry M <9 M - m) ip(x) 
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+y J(x,y)D(x - y)J(y,x) - — J 5 (x,y)D(x - y)J 5 (y,x)j. 



(2) 



This action is invariant under global chiral symmetry, i.e. ip — > e l6l5 %l) and ifj — > %l)e~ ldl5 
, provided that m = 0, and is sensible to the details of the gluon propagator D(x — y). 
Extending ()2]) to include various quark flavors is straightforward. In the following it will 
be assumed SU(2)-flavor symmetry and the current quark masses will be written as m u = 
m d = m q . 



The action ([2]) is a nonlocal version of the Nambu-Jona-Lasinio (NJL) action 25|. The 
NJL model gives a good description of low energy hadron phenomenology and it is expected 
that the nonlocal version should do a similar job, with the extra bonus of being sensible to 
the gluon propagator. The NJL model is nonrenormalizable and needs to be regularized by 
the introduction of a finite cut-off A. For the local version, where D(x — y) oc 5(x — y), 
typical cut-off values range from 600 MeV up to 800 MeV. For the model described by Eq. 
02]), besides the regularization, one also needs to know the gluon form factor D(x — y). 

First principles calculations of the gluon propagator have been performed using lattice 



QCD and DSE. The momentum space propagator 

/ o\2k— 1 



Dip 1 ) 



(p 2 + Xqcd) 



2k 
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is able to describe both the scaling (k > 0.5) and decoup 
and the lattice data up to p ~ 800Mey - see Refs. 
Aqcd stands for an infrared mass scale. 

Let us define the dimensionless form factor in momentum space as 

2 / 2 \ ^ 

/(/) = AW) = ^ L, P , 2 e{A- P ). 



(3) 



k = 0.5) infrared DSE solutions 



18[ for discussions. In Eq. 



P 2 \P 2 + A QCD 



(4) 



The constant Z in Eq. fl3]) will be included in the definition of the coupling constant G, 
which multiplies the quark currents of the nonlocal theory. G carries the dimension of a 
length squared. In f(p 2 ), A is the cut-off. In a first step we assume Aqcd = A. The form 
factor f(p 2 ), together with typical form factors considered in the literature, is shown in Fig. 

m 

The model has four parameters {G,m q , k,A} to be determined in vacuum. Given that 
Eq. ([3]) reproduces the lattice data for p < 800 MeV, we take A = 800 MeV. For each k, G 
and m q are adjusted to reproduce the experimental pion mass M n = 139.57 MeV and decay 
constant / w = 92.4 MeV. 

We also take in consideration the light quark condensate (qq). At a renormalization scale 
ofw = 1 GeV, (qq) 1 / 3 = -240 MeV, while at fi = 2 GeV it becomes (qq) 1 / 3 = -270 MeV 
27], |28(. The difference is mainly due to the running of m from 5.5 MeV, for p, = 1 GeV, to 
4.1 MeV if p = 2 GeV - see Ref. 29(. Furthermore, due to phenomenological reasons, we 
will restrict our study to constituent quark masses M q in the range 300 - 400 MeV, which 
are typical values of M q in this type of models. The decay width for n — > 77 will be used 
to distinguish the different sets of parameters. 

Our nonlocal model belongs to a class of models which have been explored by previ- 



ous authors 
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321 ] and we will not discuss how to compute the different quantities. The 



technical details can be found in [30 1 . 

We start the discussion of the results with a technical comment on the NJL-type model 
calculations. The fermionic action of NJL-type models has ultraviolet divergences and 



requires a regularizing cut-off 



33| . However, not all the integrals are divergent and some 
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4.093 
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K = 0.45 


4.146 
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3.22 


K = 0.50 


4.187 


360.5 


271.4 


6.441 


5.44 


K = 0.529 


4.205 


383.6 


271.1 


7.491 


7.79 


K = 0.548 


4.214 


400.1 


271.0 


8.263 


10.19 


K = 0.60 


4.231 


453.3 


271.4 


10.776 


26.20 



TABLE I. A collection of parameters which reproduce the experimental M v and f v for a cut-off 
A = Aqcd — 800 MeV. The results differ essentially on the value of the decay width rV-^ (last 
column). Recall that the experimental width is 7.78(56) eV. 



authors chose to regularize only the action and the divergent integrals: for the finite integrals 
they remove the cut-off. This has an important impact on r w -+ 71 . Indeed, for the local NJL 
model the width computed with all integrals re gula rized give, typically, only about 70% of 



the experimental numbers (see, for example, 3l|, |34j). In our calculation the cut-off A = 800 



MeV was used at all the stages of the computation. 

We would like to remember the reader that our work has two main goals. On one side we 
would like to define a nonlocal model for low energy physics which is, as much as possible, 
rooted in first principle computations. This motivates our choice for the form factor. On 
the other side, we would like to connect low energy physics with the behavior of the gluon 
propagator D(p 2 ) in the infrared region and, in particular, its value for p = 0. 

In what concerns M n and /„■, for a cut-off of A = 800 MeV, it turns out that the 
experimental values can be reproduced by a large set of parameters, including k above 
and below 0.5 - see Table [I] The values on Table [J show that the quark condensate is 
almost independent of k and always close to the expected value (qq) 1 ^ 3 = —270 MeV. It 
is important to point out that this model is consistent with the Gell-Mann-Oakes-Renner 
relation (GMOR) 

M 2 Jl = -2m q (qq) Q) (5) 
preserving chiral low-energy theorems and current algebra relations |26j. As matter of fact, 
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FIG. 2. yyy clS db function of k for A = Aqcd = 800 MeV. The excluded regions, k < 0.41 and 

k > 0.55, are associated with a constituent quark mass outside typical phenomenological values 
M q < 300 MeV or M q > 400 MeV. See also Table [Q 

we can calculate the GMOR value for the current quark mass 

M 2 f 2 

m GMOR = jy^K ^ (6) 



as an indicator of the validity of this low energy theorem. Taking for example K = 0.5, 
m GMOR _ 4 ^59 MeV which differs less then 1% of the calculated value m q = 4.187 MeV. 

The width r 7r _ i , 77 , with an experimental value of 7.78 (56) eV [35|, distinguishes between 
the various scenarios for D(0) - see Fig. |2j Clearly, the two pion decay points towards a 
k above 0.5 and a D(0) = 0. This conclusion is unchanged when A ^ Aqcd and for other 
cut-off values - see Fig. |3j 

The results reported in Fig. Elshow that to reproduce the experimental T^^, the cut-off 
A and the infrared mass scale Aqcd are correlated. The cut-off increases with Aqcd, with 
A and Aqcd being of the same order of magnitude. Moreover, for each Aqcd, there is only 
a narrow band of cut-off values which are able to reproduce the experimental width. For 
example, for Aqcd = 750 MeV, A has to be in the interval 761-885 MeV to reproduce 
r^-^y up to two standard deviations. Furthermore, k seems to decrease slightly when Aqcd 
increases. Given the correlation between A and Aqcd, this suggests that for higher values of 
the cut-off, in principle, one could reproduce the full set of observables for k = 0.5 or smaller. 
However, since our nonlocal model does not take into account the vector and axial-vector 
bilocal quark currents, that the light J = 1 spectrum starts at ~ 800 MeV, an upper limit of 




FIG. 3. Tn-t^ as a function of k for several values of A and Aqqd- The circles and squares mean 
M q = 300 MeV and M q = 400 MeV respectively. 

A around 800 MeV seems natural. If one wants to consider higher cut-off values, one should 
also include the vector and axial-vector currents in the model. 

The decoupling type of propagator can be investigated setting k = 0.5 in Eq. ([3]). Then 
the gluon propagator becomes 

P t m gluon 

M giuon takes the role of Aq C d an d can be interpreted as an effective gluon mass. Following 
the same procedure as before, it follows that the model for k = 0.5 also reproduces M n , f n 
and r^^, provided the gluon mass is adjusted as a function of the cut-off A. Our results 
are summarized in Table [III and Fig. HJ for various cut-off values. 

The effective gluon mass M g \ uon = Aqcd is within typical values found in the literature, 



sec 



and references there in, but shows a strong dependence with the cut-off. M g 



is a linear function of A - see Fig. HI In what concerns the quark condensate, the model 
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FIG. 4. M g i uon (A) required to reproduce the experimental T 



7T— >77 



shows that (qq) increases with M g i uon , i.e with the cut-off A. For k = 0.5, to reproduce the 
experimental value of the condensate, i.e. to have — (qq) 1 ^ 3 = 270 MeV, it turns out that 
the gluon mass is M gluon = 878 MeV for A = 800 MeV. 

Finally, we would like to refer that the GMOR relation is satisfied always, i.e. for k > 0.5 
and k = 0.5, within the same level of precision. 

In this paper, the connection between QCD and the NJL model is explored to investigate 
possible links between the quenched gluon propagator and low energy hadronic phenomenol- 
ogy. The relation between the nonlocal model and QCD used is valid at the leading order 
level. It has been shown that nonlocal covariant extensions of the NJL model have the 
advantage, over the local theory, of the next-to-leading order corrections are relatively small 
37j |. Therefore, we expect that, at least qualitatively, our results will still be valid beyond 
the leading order. 

Our main goal is to connect D(p 2 ) with pion physics. The results discussed previously 
suggest that low energy physics is essentially blind to the gluon propagator at very low 
momenta. Indeed, the results show that one is free to use either the decoupling of scaling 
scenarios for D(p 2 ), provided that the model parameters are chosen appropriately. 

For a scaling-type scenario, M n and f n are not sensitive to the deep infrared gluon 
propagator. Indeed both values can be reproduced for k above and below 0.5. However, 

yyy is sensitive to k. Curiously, for A = Aqcd = 800 MeV, the value of k which reproduces 

the experimental decay width, k = 0.529, is in perfect agreement with the infrared exponent 
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A [MeV] 


M gluon [MeV] 


m q [MeV] 


M q [MeV] 


-{qq) 1 / 3 [MeV] 


GA 2 


750 


843.8 


4.5 


460.9 


264.2 


6.03 


800 


878.7 


4.2 


409.8 


271.6 


5.39 


813 


888.5 


4.1 


400.0 


273.5 


5.27 


850 


917.7 


3.9 


373.3 


278.9 


4.96 


900 


959.1 


3.6 


345.7 


286.3 


4.66 



TABLE II. Summary of the results for k = 0.5. 



measured from 
quantization 



38 



attice QCD in 



19l . [24J and the estimate from time-independent stochastic 
. Of course, given the level of approximations used here, the good agreement 
is, probably, an happy coincidence. 

For a decoupling- type scenario, the "measured" effective gluon mass is within typical 
values. From the values of the quark condensate, the massive gluon model seems to prefer 
a M gluon ~ 878 MeV. 

The parameterizations discussed in this work gather the essential ingredients of the scaling 
and decoupling solutions of the DSE and the lattice propagators for p G [50, 800] MeV. 
Besides the link with the infrared gluon propagator, the models analyzed also provide an 
effective nonlocal chiral quarks model and, in this sense, are a good starting point towards 
the construction of a more realistic nonlocal low-energy effective theory for QCD. 
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